Stripes Induced by Orbital Ordering in Layered Manganites 
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Spin-charge-orbital ordered structures in doped layered manganites are investigated using an orbital- 
degenerate double-exchange model tightly coupled to Jahn- Teller distortions. In the ferromagnetic 
phase, unexpected diagonal stripes at x—l/m (m=integer) are observed, as in recent experiments. 
These stripes are induced by the orbital degree of freedom, which forms a staggered pattern in the 
background. A 7r-shift in the orbital order across stripes is identified, analogous to the 7r-shift in 
spin order across stripes in cuprates. At x—1/4 and 1/3, another non- magnetic phase with diagonal 
static charge stripes is stabilized at intermediate values of the t2 g -spins exchange coupling. 

PACS numbers: 75.30.Kz, 75.50.Ee, 75.10.-b, 75.15.-m 



Manganese oxides are currently attracting considerable 
attention JjJ, due to the complex interplay among spin, 
charge, and orbital degrees of freedom, which induces a 
rich phase diagram as well as Colossal Magneto- Resistant 
(CMR) properties. There are clearly two types of dom- 
inant states in these compounds. For example, in per- 
ovskite manganites such as Lai-zCa^MnOa, in the re- 
gion 0.22<a;<0.5 a ferromagnetic (FM) metallic phase is 
the ground-state at low temperature. On the other hand, 
at x>0.5, a charge-orbital-spin ordered state is stabilized. 
The competition between these two states is at the heart 
of recent theories that explain the CMR effect in man- 
ganites as arising from mixed-phase tendencies Q . 

This two-phase metal-insulator competition and con- 
comitant large MR effect occurs also at the Curie tem- 
perature Tc, at densities where a FM phase exists at low 
temperature. Several experiments have clearly revealed 
the mixed-phase characteristics of manganites near Tc 
H . While it is natural to assume that one of the com- 
peting phases is FM metallic, the properties of the in- 
sulating phase are still unknown. Recently, consider- 
able progress has been made in this context. Just above 
Tc evidence for the existence of short-range stripe-like 
charge ordering has been obtained with neutron diffrac- 
tion and X-ray scattering studies. For Lai-^CazMnOa 
j|, "diagonal" stripes, i.e. charge ordering (CO) along 
the (l,l,0)-direction g, appear for 0.2<:c<0.3, while for 
x<0.2, "bond" stripes, i.e. CO along the (1,0,0) or (0,1,0) 
direction, have been revealed. For Laa— 2xSri-|_2 X Mn207 
||, short-range bond stripes have been detected in the 
wide range 0.3<a;<0.5. These results lead to the intrigu- 
ing possibility that the insulating phase that contributes 
to the CMR near Tc may also be FM but with striped 
features, a remarkable novel result. Such a state would be 
puzzling since stripes in cuprates are associated with the 
creation of rivers of holes to avoid having the individual 
charges "fighting" against the antiferromagnetic (AFM) 
background. Thus in a FM state, where hole movement 
appears optimal, charge is naively not expected to form 



stripes, in contradiction with experiments. 

In order to understand this puzzling complex problem, 
in this Letter the two-dimensional (2D) double-exchange 
(DE) model coupled to Jahn- Teller (JT) phonons is inves- 
tigated using unbiased computational techniques. This 
model has been already successful in reproducing the A- 
type AFM state at x=0 [§, the CE-state at x=0.5 [§, 
and the complex structure of the a;>0.5 regime ||, as 
observed in experiments. However, a state as exotic as 
containing stripes in a ferromagnetic background has not 
been reported until now. The present effort focuses on 
the properties of 2D systems, since (i) studies in three- 
dimensional (3D) cases are technically far more complex, 
and (ii) stripe structures identified experimentally in bi- 
layer and 3D manganites are here observed in 2D systems 
as well. Our main result is the stabilization of FM states 
with stripe order at x=l/m (m=integer), a surprising re- 
sult whose origin lies in the concomitant orbital order. 

Regarding single-layered manganitesLai-^ri+JVInO-i, 
the undoped compound is AFM juJ, while at x~0.5 a 
CE-type AFM CO phase has been identified 0. How- 
ever, for 0.0<x<0.5, a complex "spin glass" behavior has 
been experimentally observed [Q , indicating that the 2D 
ground-state properties are basically unknown. Our re- 
sults below also indicate that stripe states with CE-like 
AFM characteristics may exist at £=1/4 and 1/3, and 
they could be important for the physics of single-layered 
compounds in non-FM regimes. 

The Hamiltonian studied here is 

H = - ^ 7 '<4 70 - rf i+a7'er ~ 2J S ' ' Sj 

ia77 / (T i 

+ Jaf 22 Si ■ Sj + A }XQupi + Q 2 iT x i + Q 3 iT z i) 

+ (l/2)^(/3Q 2 i + Q 2 i + Ql i ), (1) 

where <ii ao - (fiibrx) is the annihilation operator for an e g - 
electron with spin a in the d x i_ y i (<i 3z 2_ r 2) orbital at site 
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i, and a is the vector connecting nearest-neighbor (NN) 
sites. The first term represents the NN hopping of e g elec- 
trons with the amplitude t* t between 7- and 7'-orbitals 



along the a-direction (i* 
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V3i£ a =3f£ b =l 



for a=x and fy a = v / 3^ b = v / 3iba= 3< bb= 1 for a =y> 1 

energy units). In the second term, the Hund 
coupling Jh(>0) links e g electrons with spin Si = 

S 7 a/3 ^i-ya ""^^/ 3 ( CT are ^ ne P au li matrices) with the 
localized t2 g spin Si, assumed classical and normalized 
to I Si | = The third term is the AFM coupling Jaf be- 
tween NN t 2g spins ]l3fl . The fourth term couples e g 
electrons and Mn06 octahedra distortions, A is the di- 
mensionless coupling constant, Qn is the breathing-mode 
distortion, Q21 and Q31 are, respectively, (x 2 —y 2 )- and 
(3z 2 -r 2 )-type JT-mode distortions, p\= d i J(r di 7<T , 

T xi = EaKLx^ba + 4 h(J d ia<J ) , and T zi = £ CT (4,<x d ia<x ~ 

d\ ba diba)- The fifth term is the potential for distortions, 
where (3 is the ratio of spring constants for brcathing- 
and JT-modes, treated here adiabatically pM- 
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FIG. 1. Orbital arrangement in the CO FM phase at (a) 
1=1/2, (b) 1/3, and (c) 1/4 with A=2.0. The charge density in 
the lower-energy orbital is shown, with the size of the orbital 
being in proportion to this density, (d) Total density of states 
vs energy for x=l/3 and A=2.0. 

First let us clarify the charge and orbital structure in 
the FM phase, main issue of this paper. If lattice dis- 
tortions are not correlated in different sites, holes will be 
simply distributed as uniformly as they can to lower the 
ground state energy. Namely, the stabilization of charge 
inhomogeneous structure such as stripes in the FM phase 
requires a proper treatment of the cooperative JT effect. 
A simple way to include such an effect is to optimize the 
displacement of oxygen ions uf along the a-axis at site 
i, instead of local distortions |lq| . Results at large A for 
the FM phase are shown in Fig. 1 . Using relaxation tech- 
niques to optimize {u}'s, diagonal stripes can be clearly 
observed even in a spin FM regime, a surprising result. 
The key ingredient to understand the presence of stripes 
in a spin magnetized phase is the presence of the con- 
comitant staggered orbital order. Individual holes doped 



into the £=0 FM orbital-ordered (00) state produce a 
distortion of the ordered background. This energy lost 
is minimized if holes share the distorted regions, forming 
complex patterns such as stripes. A similar reasoning 
is usually applied to the rationalization of stripe forma- 
tion in nickelates and cuprates Jlq| , with the important 
conceptual difference that in those compounds the back- 
ground in which the stripes are created is spin AFM, i.e. 
the real spin is active, while in the present study the spin 
is FM and the orbital background is active. If the or- 
bital degree of freedom is associated to a "pseudo-spin" 
up and down, then an analogy between manganites and 
nickelates/cuprates can be established, replacing the real 
spin of the latter by the orbital degree of freedom of the 
former. In fact, with this analogy the pseudo-spin and 
charge structure at x—l/A (Fig. 1(c)) becomes the same 
as the real-spin and charge stripe structure of hole-doped 
La2Ni04, at the same hole density |l7|. In the stripe 
phase, a gap of the order of the JT energy opens around 
the Fermi level (see Fig. 1(d)). Thus, "pseudogap" fea- 
tures may be expected as precursors for stripe formation 
even at high temperatures such as Tq- 

Note that the stable charge-orbital stripes with the 
arrangement of Fig. 1 can appear when the distance d 
between diagonal hole arrays is equal to mao, where m 
is an integer and arj is the lattice constant along the a- 
axis. Since d—ao/x from Fig. 1, x should be equal to 
1/m for the appearance of stable diagonal charge stripes. 
For l/(m+l)<x<l/m, it seems possible (but at this 
early stage it is a conjecture) that a mixed phase of 
two charge-orbital ordered states with a;=l/(m+l) and 
1/m appears, consistent with the phase separation sce- 
nario H and also with recent synchrotron X-ray scat- 
tering measurements [Q. Based on this scenario, the 
charge ordering observed in experiments at a;~0.3 in the 
FM regime may be understood as a mixture of diagonal 
stripes at 2;= 1/4 and 1/3 (Figs. 1(b) and (c)), if those 
patterns are assumed to be stacked along the z-axis due 
to the influence of Jaf, a reasonable assumption based 
on previous x=0.5 calculations B. In the orbital correla- 
tion function T z (q) =Xa,j e * q ( r zi r zj) ; peaks appear 
at (7T±<5 m , 7r±f5 m ) with 5 m =(l— 2/m)7T. The deviation 
from (7r,7r) for m>2 is caused by a tt- shift in the orbital 
order across the stripe (Fig. 1), and it can be informally 
referred to as "orbital incommensurability" by analogy to 
the spin incommensurability found in cuprates and nick- 
elates. It is important to remark that the present ideal- 
ized charge stripes will likely be destabilized by thermal 
and/or quantum fluctuations. Thus, in actual materi- 
als, it is expected that the stripes will be dynamical as 
in cuprates and only vestiges of stripes may be detected 
p9[ , together with pseudogap features, consistent with 
the phase separation tendency for l/(m+l)<x<l/m. 

Now let us consider the charge-orbital structure in the 
AFM phase by using non-cooperative JT phonons. This 
will allow us to report yet another striped state which 
could be observed experimentally, this time with an over- 
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all zero magnetization. As found at x=l/2 @], the AFM 
phase such as CE-type state is not sensitive to the JT 
phonons treatment, since there exists a strong constraint 
due to the DE mechanism for the e g electron kinetic mo- 
tion, masking differences in the character of JT phonons. 
Note, however, that both the local lattice distortion and 
t2g spin direction should be determined independently at 
each site by optimizing the total energy. Using relaxation 
techniques to optimize {Q}'s and {S}'s at fixed electronic 
density p0[) , the phase diagram at x—1/4 has been here 
obtained (Fig. 2(a)). Its overall features can be under- 
stood from the competition between the e g electron ki- 
netic energy and magnetic energy of t2 g spins regulated 
by Jaf- At small Jaf the system becomes FM to im- 
prove hole movement, and a metal-insulator transition 
occurs at a critical value of A, separating FM CO and 
FM charge-disordered (CD) states. In the other limit of 
large Jaf, an AFM phase is stabilized since the mag- 
netic energy among t 2g spins dominates. The most in- 
teresting result of Fig. 2(a) is that at intermediate val- 
ues of Jaf, a novel spin-ordered state analogous to the 
CE-type phase at £=1/2 has been found. The complex 
optimized spin pattern is shown in Fig. 2(b). A simi- 
lar CE-like spin arrangement is also found at x=l/3 (see 
Fig. 2(c)) pljj . These configurations are here called the 
"zigzag" AFM (Z-AFM) states, since t 2g spins form one- 
dimensional (ID) zigzag paths where e g kinetic energy is 
gained, stacking with anti-parallel spins in the direction 
perpendicular to those paths to gain magnetic energy. 
Since this Z-AFM phase can take partial advantage of 
both energies, it is reasonable that it is stabilized in be- 
tween the FM and AFM phases. 
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FIG. 2. (a) Zero-temperature phase diagram at £=1/4 and 
large Jh (see text for notation). The metallic FM-CD and 
insulating FM-CO boundary is obtained monitoring the aver- 
age JT distortions with increasing A. (b) Ground-state spin 
pattern on an 8x8 cluster for intermediate values of Jaf at 
a:=l/4 (Z-AFM phase in (a)). Open and solid arrows indicate 
up and down spins, respectively, (c) Spin arrangement on the 
6x6 cluster in the Z-AFM regime at £=1/3. 

Consider now the charge and orbital structures of the 
phases in Fig. 2(a). In the FM phase, for the values of 
A investigated, an OO phase appears (not shown), irre- 



spective of the charge structure. As mentioned above, for 
non-cooperative JT phonons, holes are distributed uni- 
formly to lower the ground state energy, compatible with 
the concomitant orbital order. In the Z-AFM phase, the 
charge-orbital arrangements for large A are schematically 
shown in Figs. 3(a) and (b) for x=l/4 and 1/3, respec- 
tively. The local charge density is found to be constant 
along the diagonal directions denoted by the broken lines 
in Fig. 3(a). This result indicates a tendency toward the 
formation of diagonal charge stripes in the Z-AFM phase 
for 2:= 1/3 and 1/4. In Fig. 3(c), the local charge den- 
sities on the diagonal lines, n(S\), 11(62) , and ^(#3), are 
shown vs A. At small A, diagonal charge stripes are ob- 
served, concomitant with a peak around q=(7r/2,7r/2) in 
the charge correlation function n(q) p2]| . On the other 
hand, at large A, 71(83) is very small, while n(8i) and 
n(S 2 ) are almost unity. Namely, holes are mainly located 
along the line £3 , indicating the formation of a clear diag- 
onal charge stripe pattern (see Fig. 3(a)), and the inten- 
sity of the (7r/2, 7r/2) peak in n(q) becomes larger than 
at small A. Note that in the Z-AFM phase the e g elec- 
tron motion is restricted to the ID zigzag FM chains due 
to large Jh- When a finite electron-lattice coupling is 
introduced in this ID system, a Peierls instability should 
occur, leading to a charge-density-wave state. Since in 
the Z-AFM phase the same chain is simply stacked along 
the diagonal direction in the x-y plane, diagonal charge 
stripes occurs naturally [ p3| . 

To understand the shape of the zigzag chain (Fig. 2(b)), 
consider the limit of A=0 ||||. Even without phonons, 
by straightforward diagonalization it can be shown that 
the zigzag chains have a spectra corresponding to a band- 
insulator due to the periodic changes in hopping ampli- 
tudes along zigzag chains (t^ h — — i* h ), which induce gaps 
in the energy spectra ||. At 2;= 1/4, there are nine inde- 
pendent possible types of zigzag chains on an 8x8 lattice 
with periodic boundary conditions. Let us label these 
zigzag chains using "bits", and 1, representing the x- 
and y-directions, respectively. The nine possible config- 
urations are in Fig. 3(d). The zigzag chain of Fig. 2(b) 
is given by the periodic sequence "00101101", and its en- 
ergy is the lowest among the possible candidates, since it 
provides the largest bandgap. When the electron-phonon 
interaction is adiabatically switched on, it is reasonable 
that the AFM phase Fig. 2(b) would still be stable for 
a finite A range pi} ]. Note that in Fig. 3(d), excited 
states characterized by other zigzag paths have small 
excitation energies, such as 0.001~0.01 in units of £* a . 
Thus, at temperatures as low as a few degrees Kelvin the 
ground-state Fig. 2(b) can be easily distorted into other 
zigzag spin patterns. This fragility of the ground-state 
may be related with the "spin glass" features observed in 
single-layer experiments |l2| , leading to an overall orbital 
disordered state. However, note that seven of the nine 
competing states, as well as the Z-AFM phase, have diag- 
onal stripes. Then, even in a mixture of these states, the 
stripe direction is not random. Thus, it is possible that 
indications of diagonal charge stripes may be present in 
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the "spin glass" phase of single-layer manganites. 
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FIG. 3. (a) Schematic view of spin-charge-orbital order- 
ing in the Z-AFM phase at large A and £=1/4. Open and 
solid symbols indicate up and down t2 S spins, respectively. 
The lobes indicate d 3x 2_ r 2- or d 3H 2_ r 2-orbitals in the Mn 3+ 
site, and the small circles denote the hole sites in which 
d a .2_ !/ 2-orbital is occupied by a tiny amount of charge, (b) 
Same as (a), but for £=1/3. (c) Charge densities (see (a)) vs 
A at a;=l/4 in the Z-AFM, where results at a given A are Jaf 
independent, (d) Energies of several paths for £=1/4, A=0, 
and Jaf =0.25 in the Z-AFM regime (see text for notation). 

In summary, novel striped charge-orbital ordering has 
been found in realistic models for manganites. Diagonal 
stripes in the FM phase have been observed at densi- 
ties x=l/m, with m an integer, and also in the CE-like 
Z-AFM phase. The orbital degree of freedom orders in 
a x—0 staggered pattern in between the stripes, play- 
ing a key role in stabilizing these structures, similarly as 
the real spin does for stripes in nickelates and cuprates. 
Our results have implications for the recently discovered 
short-range charge ordering effects in neutron scattering 
experiments, as well as for future experiments, particu- 
larly involving layered manganites. 
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[1] See, e.g., Colossal Magnetoresistive Oxides, edited by Y. [22] 
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Once the stripes melt, together with the dynamical as- 
pect of the charge, the orbital order will no longer be 
long-ranged but short. Then, melted stripes may lead to 
orbital liquid properties (M. Belesi et al, preprint). 
To reduce the effort in our calculations, still capturing the 
main effect of JT distortions, /3 is set to oo effectively sup- 
pressing the breathing mode, while Jh=20 (but results 
at Jh=oo are essentially the same). 

At x=l/3 and large A, the zigzag FM chain (Fig. 2(c)) 
is divided into small FM clusters for non-cooperative JT 
phonons, while for cooperative case the Z-AFM phase 
(Fig. 2(c)) is stable. For large A, the cooperative effect 
is crucial to stabilize the Z-AFM structure, as well as to 
observe the diagonal charge stripes in the FM phase. 
Note that the charge density is nearly uniform at inter- 
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mediate A such as 1.2 in Fig. 3(c). In this interesting 
regime, a peak in the orbital correlations exist, while the 
analog in the charge sector is not prominent. Such an OO 
but CD £=1/4 state may be related to those observed in 
recent exper iments for Pri-^Ca ^MnOs (M. v. Zimmer- 
mann et ai, 2ond-mat/0007231). 

[23] In the Z-AFM phase, the charge pattern depends on A 
and the shape of the FM path, but the orbital arrange- 
ment is mainly determined by the zigzag geometry. 

[24] Note that the results of Figs. 2(b,c) and the analysis of 
competing zigzag paths have been obtained assuming ei- 
ther FM or AFM links among localized spins, reasonable 
at large A. If the *2 g spin directions are optimized for 
small A, spin patterns similar to Fig. 2(b) and (c) are 
stabilized, but with spin directions slightly disordered. 
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